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Abstract 

PU . Physics beyond the Standard Model could be measured indirectly, through its 

effects on Standard Model observables. One place to look for such effects is in 
the semileptonic decays of B mesons. In order to constrain the possible role of 
new physics on semileptonic B decays, we present the most general low energy 
effective Lagrangian constructed from all dimension six four-fermion operators 
that contribute to the process b — > due. We then use it to compute the 
^^ I corrections due to new physics to the differential decay rates for the exclusive 

Q I processes B — > D^*>£i?£, as well as the inclusive decays B — > Xcii'i. Both 

0^ ' inclusive and exclusive rates are expressed in terms of a set of parameters 

that characterize the types of four-fermion interactions that can be induced 



by physics beyond the Standard Model. Although it is not particularly useful 
to carry out a full analysis until data from the next generation of B factories 



Q I becomes available, here we illustrate how the existing experimental results 

i-^ ' can be used to constrain some of these parameters. 
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I. INTRODUCTION 

Despite the fact that the Standard Model successfully accounts for most of the observed 
experimental results (see [|l| for a review), it is incomplete, and must be replaced by a 
more fundamental theory at some high energy scale. Whatever form physics at a higher 
scale might take, it will generate low energy effective couplings that contribute to Standard 
Model processes. These couplings provide the opportunity to search for signatures of new 
physics indirectly, by measurements of the deviation of experimental observables from their 
Standard Model predictions. 

Rather than compute the effect of physics beyond the Standard Model within the frame- 
work of a specific model, our ignorance dictates that we include all effective interactions 
consistent with the symmetries of the Standard Model. The resulting observables are ex- 
pressed in terms of the parameters of the generalized interaction, which can be constrained 
by the experimental data. Since these parameters characterize the most general interaction, 
their experimental bounds can be used to rule out or guide in the construction of specific 
models. 

Anticipating the precision data from the upcoming generation of B factories, one place to 
look for the effects of new physics is in observables related to B meson physics. Therefore, in 
this paper we set up a parametrization for a set of observables associated with semileptonic B 
decays into charmed mesons. Specifically, in Section |I| we construct the low energy effective 
Lagrangian for the process b -^ ciui induced by new physics. We do so by including all the 
dimension six, four-fermion Lorentz scalars that contribute to the process, but excluding 
higher dimension operators. We also introduce a set of real parameters bilinear in the 
coupling constants of this effective Lagrangian. These parameters can be used to characterize 
the low energy effects of the new physics. 

Using the generalized interaction of Section ||, we then calculate at leading order in 
HQET 0-§| the contribution from new physics to the differential decay rates for B — >■ 
Diui and B -^ D*iui. The results for these observables, expressed in the parametrization 
of Section || appear in Section |T| and Section |^ respectively. The parametrization of 



Section |I| also appears when we consider the corrections induced by physics beyond the 
Standard Model on the inclusive semileptonic decay of B mesons into final states containing 
a c quark, B -^ Xcdvi. In Section |V|, these deviations are calculated at leading order in 
perturbation theory, which is justified for m^ ^ Aqcd 0- In Section |V| we use the results 



of the previous sections to discuss briefly the constraints which the existing experimental 
data imposes on the parametrization of Section y. 

Finally, we note that new physics effects on observables related to 6 -^ c^vi have been 
studied extensively for restricted types of non-Standard model interactions. The decays 
B -^ Xiui (specially for decays into r-leptons) have been used to constrain classes of models 
with scalar interactions pj-p!^, such as those which arise in models with extended Higgs 
sectors. Other constraints on such types of models have also come from observables related 
to the exclusive decays B -^ Diui [1T3|-[T5[] as well as B ^ D*iu£ [^,|lj,|l^. The process 
b -^ clvi has also been used to analyze possible vector and axial vector couplings beyond the 
Standard Model V — A. In particular, model-independent analyses of right-handed b —>■ c 



quark current contributions to 5 — > Xiui and B -^ D*ivi were performed in |T^,|T8[. These 



were motivated by an extension of the Standard Model, the Left-Right Symmetric model [[19 



which contains gauge bosons that couple to right-handed fermions. Also inspired by this 



class of models were ||T3| , |20[1 , who constructed observables based on the exclusive decays 
B — > D^*Hu£ in order to constrain right-handed vector couplings. The combined effects of 
non-Standard Model scalar, pseudoscalar, vector, and axial vector interactions were studied 
in [^ within the context of CP violating polarization observables related to B ^ Dtv^ 
and B -^ D*iui. 

II. DEFINITIONS 

Our starting point for the calculation of semileptonic B decays to charmed mesons is the 
following interaction Hamiltonian: 



'Lint — 't-SM ~r 



4Gf 



V2 •™^-'^^ 



K6E<e[cr^&/.][4rv,]+h.c. 
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TisM is the effective V — A Hamiltonian which mediates the Standard Model Weak inter- 
action process b — > ciui for energy scales much less than mw- The second term represents 
the low-energy effective interaction generated by new physics at some energy higher than 
rriw Rather than adopt a particular model for this new physics, we include all the lowest 
dimension operators which contribute to 6 ^ ciui at tree level and which also respect the 
Lorentz and gauge invariance of the Standard Model. These are the dimension six four-fermi 
operators which appear in Eq. (|I]). One way to organize their contribution \T2] is to sum 
the index 7 over S, V, T, with 



r^ = 1; (r^)M = Y; {r^r = ^'^ = ^[7^71, (2) 

while projecting the 6-quark and lepton I into fields of definite chirality (respectively denoted 
by h^ and f^, with /i, e G {L, R}) and summing over all values of /i, e. Given 7, the chirality 
of the c quark is fixed by /i, and the chirality of the neutrino field vi is fixed by the value of 
e. Note in particular that we are including operators which contain a right-handed neutrino 
field. These occcur for 7 = S', T and e = L, or 7 = V^ and e = R. 

The complex coefficients f?^^ are a set of twelve dimensionless coupling constants (scaled 
by the Fermi constant Gp), one for each choice of 7, /i, and e. (In the abscence of right- 
handed neutrinos, g^]^ = g^j^ = gj^j^ = 0, and the number of coefficients would be reduced 
to six.) If written in terms of the fields that transform as definite representations of the 
Standard Model SU{3) x SU{2) x f/(l) gauge group, no insertions of the Standard Model 
Higgs doublet are needed to make the dimension six operators of Eq. (|l]) gauge invariant (we 
take the right-handed neutrino to be a Standard Model gauge singlet). Dimensional analysis 
then implies that gj^^ ~ (9(?7i^/M^), where M is the energy scale of the new physics. 

Only certain real combinations of these coupling constants will appear in the expressions 
for the calculated observables. To first order in gj^^, we have 

l3o = 2Re{glj], 
P'^ = 2Re{gl^], 
eo = 2Re{(7f^}, 
e[, = 2Re{^|^}, 
Po = 2Re{glj,}, (3) 

which in principle could be present even in the abscence of a right-handed neutrino. However, 
the last three of these can only appear for non-zero lepton mass. We expect these five 
parameters, which arise by interference with Tismi to give the most significant contribution 
to signals of new physics. Eq. (^ includes only the interference with the leading order part 
of the Standard Model term. In reality, it should be modified by higher order Standard 
Model effects (such as radiative corrections), which we ignore here. 

The remaining parameters are second order in the gj^^ and are suppressed in comparison 
to the previous set. We classify them by the following types: 

• Scalar-scalar: 
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Vector-vector: 
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P" = 2Re{gl,gl,* + gl^gl^*} 



(5) 



Tensor-tensor: 
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Scalar-tensor 



5 = 2Re[glL9RL + 91r91r} , 
5' = 2Re[glLglj^* + gl^gl,^*]. 



(7) 



Scalar- vector (only appear for non-zero lepton mass): 

e = 2Re {g^LR* + 91l9rr* + 91r91l + ^L^Il*} 
e' = 2Re {^^^f^* + ^L*?!^* + ^L^L* + ^L^Il*} 



(8) 



Tensor- vector (only appear for non-zero lepton mass): 



^LL^Li? +9rr9rL 

p = 2Re [gRL9LR + 9lr9rl 






(9) 



III. B -^ Dluf 



As in [^, we write, in the rest frame of the i?-meson: 
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1 + cos 0)2 \H+r + (1 - cos 0)2 \H_r + 2sin2 9 \Hn 



(10) 



where pn is the D three-inoinentum in the B rest frame, 6 is the angle between the lepton and 
the D meson in the rest frame of the 'virtual PF-boson,' (i.e, the frame in which pe = pn) 



and g^ = {pb —pnY- \H±\ and \Hq\ are helicity amplitudes, which we decompose as: 
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In these expressions, H^'^^ and Hq^' are the helicity amplitudes computed assuming 
only the Standard Model V — A coupling. The terms Aj_ , Ag are contributions from 
the interference of the Standard Model part of the amplitude (taken only to leading order) 
and other terms in Eq. (p. Finally, the terms Aj_ and Ag are second order in the new 
interactions of Eq. (|l]). Here we calculate the interference and the second order terms by 
matching hadronic matrix elements of the quark-quark operators in Eq. (|l]) to HQET matrix 
elements, neglecting both perturbative 0{as) and non-perturbative OII/itlq) corrections. 
Higher order corrections to the Standard Model contribution to the decays B — > D^*Hh'(, can 
be found in [p^ . 

Writing w = {m^ + m]^ — q^)/2mBmD, and denoting the Isgur-Wise function p by ^(w), 
we find for the interference tems: 
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and 
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(13) 



Likewise, the second order terms are 
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(14) 



(15) 



where A^ and A^ represent corrections due to non-zero lepton mass (which are only non- 
neghgible for decays into r-leptons): 
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IV. B 



D* 



As in the previous case, we write in the B rest frame 
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The kinematic variables Pd*? & and q^ are defined as before (just replace D with D* in the 
definitions of the previous section). Also as before, we separate the helicity amplitudes as 
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Once again, the interference terms only take into account the leading order part of the 
Standard Model amplitude. Calculating these as well as the second order terms by matching 
hadronic matrix elements onto HQET ones (neglecting 0{as) or 0{l/mQ) corrections), we 
find: 
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The second order terms are 
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with the following corrections due to the effects of a non-zero lepton mass: 
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V. B ^ Xlut 

Starting from the operator product expansion (OPE), it can be shown that in the limit 
rrib ^ Aqcd the inclusive semileptonic B decay rate is equivalent to the perturbative quark 
level h decay rate . This makes it possible to calculate the effects of new physics on the 
inclusive decay B —>■ Xcii^e directly from Eq. (|l]), by calculating the rate for the quark level 
process b -^ ciui. We write the differential decay rate as 
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(SM) \~-i ^ "/ (1) \""i ^^t^^v / (2) 

where g^ = [pb—PcY and Ei, Ey are the h quark rest frame energies of the lepton and neutrino 
respectively. Eq. (p6D is not meant to imply that the differential rate is a contribution from 
three distinct processes. Instead, we are merely decomposing the observable into a piece 
arising solely from Standard Model physics, an interference term, and a term that is second 
order in the new physics. At present, nonperturbative corrections to the Standard Model 
term are known to order (AqcD/TTife)^ [^], while only the part proportional to a^/3o of the 



NLO perturbative corrections is known ||25| (/3o being the coefficient of the one-loop QCD 
beta function). Here, we calculate the last two terms in Eq. (p6D, neglecting 0{as) or 
O^K^Qj^/mk) corrections. Introducing the kinematic variables 
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and y = Eg/Eg^^^ with Ee 



{ml — ml + m,l)/2m,b, as well as x = m,e/m,b and z = m.c/mb, 



we find, in the rest frame of the b quark 
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which, as in the exclusive decays, includes only the interference with the leading order 
Standard Model contribution. The remaining contributions are of second order in the new 
physics: 
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Integration of these expressions over E^ is trivial. Performing the q^ integral over its physi- 
cally allowed region gives the lepton spectrum 
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is the g^ integral of Eq. (|28D and 
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is that of Eq. (|29|) . In these last two expressions, the coefficient functions are defined by: 
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VI. DISCUSSION 

Any effects of pliysics beyond tlie Standard Model sliould become more apparent witli 
the improved data from tlie next generation B factories, some of which are scheduled to 
go on line in the near future. Therefore, a detailed extraction of our parametrization of 
new physics should await these results. For completeness, however, we will use the formulas 
derived in the previous sections together with the existing data to obtain bounds on some 
of the parameters of Section ||. Our analysis is only meant to be illustrative and therefore 
should not be taken too seriously. 

One place to look for constraints on the parametrization of Section y is in the extraction 
of \Vcb\ from the exclusive decays B -^ D^*^ei'e- Because in the zero recoil limit {w — > 1), 
the 0{l/mQ) HQET corrections to B ^ D'^*Hh'i vanish [^, and because ^(1) = 1 0, 



these observables provide a theoretically clean way of extracting the CKM matrix element 
\Vcb\ from experiment p7| , |25|] . In the abscence of new physics, we expect the value of \Vcb\ 
obtained from B —>■ Dev^. to agree with that from the decay B —>■ D*eUe- Therefore, a 
discrepancy between the two measurements of \Vcb\ could be used to put constraints on the 
parametrization of Section |I|. 
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Integrating Eq. (0) and Eq. ([T8| ) over cos^, and making a change of variables from q^ 
to w, the exclusive decay rates can be written in the form [Q 
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where r = mr)/mB and r* = mr)* /mB- From the results of Section |T| and Section |I^, we 
find 
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Note that although the decay rate for B -^ Deve is well-behaved at w = 1, the quantity 
J-'oiw) actually has a pole there. This is only because the scalar contribution to the decay 
rate (the terms involving a and a') vanishes more slowly than [w"^ — 1)^/^ as w ^ 1. In fact, 
the observation of this behavior as w ^ 1 in the experimental data for J^£,{w) could be seen 
as evidence for the possibility of scalar contributions to the process b -^ ciue. 

We expect the second order parameters in Eq. ( pSD and Eq. ( |39D to be suppressed by a 
factor of {rriyir/My in comparison to the interference terms. Ignoring their contribution, we 
find in the zero recoil limit: 



•^D-(l)^ 
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D 
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The CLEO Collaboration has used data on the exclusive decay B -^ Dev^ at zero recoil to 
extract |K;b|J^i?(l) = 0.0337 ± 0.044 ± 0.048J:[!;[![]i |3. From the data on 5 ^ D*eve, it 
has also found |K6|^d-(1) = 0.0351 ± 0.0019 ±0.0018 ± 0.0008 [|8[. If we neglect 0(l/m^) 
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corrections, J-'j^ {1) = JF^* (1) = 1 (the l/mg corrections vanish automatically by Luke's 



theorem [26]). Eq. (^) then reads 

0.96 ^ : ' 7 ' "^n (41) 






which to lowest order in the parameters, gives (3^ ~ —0.02. 

We can also use the data on the inclusive decays to put constraints on the parameters 
of Section ^ The OPE for the lepton spectrum does not agree locally with the physical 
spectrum near the maximum lepton energy. Therefore, the result of Section |V| can only be 
compared with experiment by constructing an observable which integrates the OPE result 
over a sufficiently large region. A suitable observable (introduced in [^) for comparison of 
the OPE with the inclusive data is given by 

_ Ii.5Gey Ee^/Ee 
^ ~ ~1 rfr .p > (42) 

Jl.SGeV dEt^^ 

which is independent of the value of | Vcfe|. Using the inclusive B -^ Xeve data from CLEO, 



the authors of E^ extracted a central value of -Ri = 1.7831 GeV. (This includes contributions 



from h -^ ueVei which introduce an error of only a few percent.) If we split -Ri into a piece 
coming from Standard Model physics alone, and a correction 5Ri due to new physics, we 
find 

.„ Ii.5GeY Ee-§-dEe Rf'^'' J^^^^y -§-dEi 

'^^1 - Z /.P\ ,^ - /.rx ,^ > (43) 



/l.SGeV (^)3^ dEe /i.sGeV {^) ,^ ^E, 



,(SM) 



where Ri is the Standard Model contribution to Ri. For the purposes of this crude anal- 
ysis, we can replace it by the full Ri with negligible error. Similar reasoning allows us to 
approximate the Standard Model part of the lepton spectrum by its leading order contribu- 
tion. Using Eq. (|3l|) and Eq. (H) to do the integrals of Eq. (|4|), 



8Rx ~ -0.001a - 0.002a' - 0.005(/5o + /?) - 0.036/3' 

+ 0.004(/5^' + 13") - 0.031(5 - 0.3097, (44) 

where we have used the values irth = 4.8 GeV and rric = 1.4 GeV [^. Note, however, that 
the magnitude of the coefficients in Eq. d^ ) is rather sensitive to the particular values of rrib 
and z = mc/rrii,. Therefore, the bounds on new physics that are obtained from this equation 
will be highly dependent on the numerical values of the b and c quark masses that are used 
to evaluate them. 
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In |^9|, the Standard Model 0{l/ml) corrections to the theoretical value of i?i are 



estimated to be 

SR^ = -(0a¥ + 5.5AAi + 6.8AA2 + 7.7pi)/m|, (45) 

where ttib = {mB + Stub*)/ 4: is the spin-averaged B meson mass, which we take to be 
rfiB = 5.31 GeV [^]. Assume that the first three terms of this equation are fixed by the 
values A = 0.39 GeV, Ai = -0.19 GeV^ A2 = 0.12 GeV^ of @. The term involving pi has 



larger uncertainty. Varying it between pi = and the estimated value pi ^ (300 MeV)^ 
of [^ gives 

- 8.5 X 10-^ GeV < 6R1 < 5.3 x 10"^ GeV. (46) 

Without a better knowledge of pi , a contribution to -Ri from new physics in the range of 
Eq. (|4^) cannot be excluded, even if the other HQET parameters were known to arbitrary 
precision. Therefore, assuming that the corrections to the theoretical value of -Ri from new 
physics also lie in this range, we can derive constraints on our parameters by comparing 
Eq. ( p6D with Eq. (^). We will do so by choosing only one parameter from Eq. (^) to 
be non-zero at a time. Then the bounds are derived by taking the ratio of Eq. ( P6| ) to 
the coefficient of that parameter in Eq. (|3). Of course, this is not entirely satisfactory. 



since some of these parameters are not independent of each other. For instance, setting 
a = forces a' to vanish as well. However, this procedure should be sufficient if all we are 
interested in is a rough numerical estimate. Thus, given the uncertainty on pi, we cannot 
rule out the following range for the scalar terms: 

< a < 0.72; - 0.26 < a' < 0.41, (47) 

Likewise, for the vector parameters 

- 0.11 </3o< 0.18; |/3|<|/3o|, 

0<P'< 0.023, 
-0.20 </3j< 0.13; |/5"|< 1/3^1, (48) 

and finally —0.017 < S < 0.027 and < 7 < 0.003. Note in particular that the bounds on 
/3q are consistent with the value derived from the exclusive data. 
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